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Summary. Analogue spacetimes are powerful models for probing the fundamen- 
tal physical aspects of geometry — while one is most typically interested in ulti- 
mately reproducing the pseudo-Riemannian geometries of interest in general rela- 
tivity and cosmology, analogue models can also provide useful physical probes of 
more general geometries such as pseudo-Finsler spacetimes. In this chapter we shall 
see how a 2-component Bose-Einstein condensate can be used to model a spe- 
cific class of pseudo-Finsler geometries, and after suitable tuning of parameters, 
both bi-metric pseudo-Riemannian geometries and standard single metric pseudo- 
Riemannian geometries, while independently allowing the quasi-particle excitations 
to exhibit a "mass" . Furthermore, when extrapolated to extremely high energy the 
quasi-particles eventually leave the phononic regime and begin to act like free bosons. 
Thus this analogue spacetime exhibits an analogue of the "Lorentz violation" that 
is now commonly believed to occur at or near the Planck scale defined by the in- 
terplay between quantum physics and gravitational physics. In the 2-component 
Bose-Einstein analogue spacetime we will show that the mass generating mechanism 
for the quasi-particles is related to the size of the Lorentz violations. This relates 
the "mass hierarchy" to the so-called "naturalness problem" . In short the analogue 
spacetime based on 2-component Bose-Einstein condensates exhibits a very rich 
mathematical and physical structure that can be used to investigate many issues of 
interest to the high-energy physics, cosmology, and general relativity communities. 



2 Silke Weinfurtner, Stefano Liberati, and Matt Visser 

1 Introduction and motivation 

Analogue models of curved spacetime are interesting for a number of rea- 
sons [1]: Sometimes the analogue spacetime helps us understand an aspect 
of general relativity, sometimes general relativity helps us understand the 
physics of the analogue spacetime, and sometimes we encounter somewhat un- 
usual mathematical structures not normally part of the physics mainstream, 
with the payoff that one might now develop new opportunities for exploiting 
the traditional cross-fertilization between theoretical physics and mathemat- 
ics [2, 3, 4]. 

Specifically, in this chapter we will discuss an analogue spacetime based 
on the propagation of excitations in a 2-component Bose-Einstein conden- 
sate (BEC) [5, 6, 7, 8, 9, 10]. This analogue spacetime has a very rich and 
complex structure. In certain portions of parameter space the most natural 
interpretation of the geometry is in terms of a specific class of pseudo-Finslcr 
spacctimes, and indeed we will see how more generally it is possible to asso- 
ciate a pseudo-Finsler spacetime with the leading symbol of a wide class of 
hyperbolic partial differential equations. In other parts of parameter space, 
the most natural interpretation of the geometry is in terms of a bi-metric 
spacetime, where one has a manifold that is simultaneously equipped with 
two distinct pseudo-Riemannian metric tensors. Further specialization in pa- 
rameter space leads to a region where a single pseudo-Riemannian metric 
tensor is encountered — this mono-metric regime corresponds to Lorentzian 
spacetimes of the type encountered in standard general relativity and cos- 
mology [11, 12, 13, 14]. Thus the analogue spacetime based on 2-componcnt 
BECs provides models not just for standard general relativistic spacetimes, 
but also for the more general bi-mctric , and even more general pseudo-Finslcr 
spacetimes. 

Additionally, the 2-BEC system permits us to provide a mass-generating 
mechanism for the quasi-particle excitations [5, 6] . The specific mass-generating 
mechanism arising herein is rather different from the Higgs mechanism of the 
standard model of particle physics, and provides an interesting counterpoint 
to the more usual ways that mass-generation is achieved. Furthermore, at 
short distances, where the "quantum pressure" term can no longer be ne- 
glected, then even in the mono-metric regime one begins to see deviations 
from "Lorentz invariance" — and these deviations are qualitatively of the 
type encountered in "quantum gravity phenomenology" , with the interest- 
ing property that the Lorentz violating physics is naturally suppressed by 
powers of the quasi-particle mass divided by the mass of the fundamental 
bosons that form the condensate [7, 8, 9, 10]. So in these analogue systems 
the mass-generating mechanism is related to the "hierarchy problem" and 
the suppression of Lorentz-violating physics. The 2-BEC model also allows 
us to probe the "universality" (or lack thereof) in the Lorentz violating sec- 
tor [7, 8, 9, 10]. More generally, as one moves beyond the hydrodynamic limit 
in generic pseudo-Finsler parts of parameter space, one can begin to see hints 
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of geometrical structure even more general than the pseudo-Finsler geome- 
tries. 

While we do not wish to claim that the 2-BEC analogue spacetime of this 
chapter is necessarily a good model for the real physical spacetime arising 
from the putative theory of "quantum gravity" (be it string-model, loop- 
variable, or lattice based), it is clear that the 2-BEC analogue spacetime is an 
extraordinarily rich mathematical and physical structure that provides many 
interesting hints regarding the sort of kinematics and dynamics that one might 
encounter in a wide class of models for "quantum gravity phenomenology" . 
This is the fundamental reason for our interest in this model, and we hope we 
can likewise interest the reader in this system and its relatives. 



2 Theory of the 2-component BEC 

The basis for our analogue model is an ultra-cold dilute atomic gas of N 
bosons, which exist in two single-particle states |A) and |B). For example, we 
consider two different hyperfine states, \F = l^mp = —1) and \F = 2, nip = 
2) of 87 Rb [15, 16]. They have different total angular momenta F and therefore 
slightly different energies. That permits us, from a theoretical point of view, 
to keep mA ^ ^b, even if they are very nearly equal (to about one part in 
10 16 ). At the assumed ultra-cold temperatures and low densities the atoms 
interact only via low-energy collisions, and the 2-body atomic potential can 
be replaced by a contact potential. That leaves us with with three atom- 
atom coupling constants, C/aa, Ubb, and Uab, for the interactions within and 
between the two hyperfine states. For our purposes it is essential to include 
an additional laser field, that drives transition between the two single-particle 
states. 

In Fig. 1 the energy levels for different hyperfine states of 87 Rb, and pos- 
sible transitions involving three-level processes are schematically explained. 
A more detailed description on how to set up an external field driving the 
required transitions can be found in [17]. 



2.1 Gross Pitaevskii equation 

The rotating-frame Hamiltonian for our closed 2-component system is given 
by: 4 



4 In general, it is possible that the collisions drive coupling to other hyperfine 
states. Strictly speaking the system is not closed, but it is legitimate to neglect 
this effect [18]. 
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Fig. 1. The horizontal lines indicate the hyperfine states of 87 Rb. The arrows rep- 
resent two laser fields — with the two frequencies fi\ and Q% — necessary to drive 
transitions between the two trapped states \F = 1, vtif = — 1) and \F = 2, uif = 2), 
where the frequency difference corresponds to the energy difference of the two hy- 
perfine states. This is realized by a three-level atomic system, because the hyperfine 
states must be coupled over an intermediate level, that has to lie somewhat below 
the excited |e) states, as indicated by A. 



with the transition rate A between the two hyperfine states. Here <Ji(r) and 
4?l (r) are the usual boson field annihilation and creation operators for a single- 
particle state at position r, and a x is the usual Pauli matrix. For temperatures 
at or below the critical BEC temperature, almost all atoms occupy the spatial 
modes "^(r) and \P#(t). The mean-field description for these modes, 

ft 2 

-V 2 + Vi- ih + Uu \&i\ 2 + Un m 



2 m. 



iPi + AJP,-, (2) 



are a pair of coupled Gross-Pitaevskii equations (GPE): — > (A, B) 
(M)-(B,A). 



2.2 Dynamics 

In order to use the above 2-component BEC as an analogue model, we have 
to investigate small perturbations (sound waves) in the condensate cloud. 

5 The perturbations amplitude have to be small compared to the overall size of the 
condensate could, so that the system remains in equilibrium. 



Analogue spacetime based on 2-component Bose-Einstein condensates 5 

The excitation spectrum is obtained by linearizing around some background 
densities pio and phases iO , using: 

$i = Vp l o+ep. ll e^« +£9 ^ for i = A, B . (3) 

To keep the analysis as general as possible, we allow the two initial background 
phases to be independent from each other, and define 

S AB = 6» A0 - 0bo, (4) 

as their difference. 

A tedious calculation [5, 6, 7] shows that it is convenient to introduce the 
following 2x2 matrices: An effective coupling matrix, 

£=~+X, (5) 

where we introduced the energy-independent matrix 



This matrix contains the quantities 



Uaa Uab 
Uab Ubb 



(6) 



A cos c$ab V^aoPbo 1 
uaa = uaa ^ — , {<) 



2 Pao 
Acos(5 A b v'paqPbo 1 

"2 'pJo' 



Ubb ^ Ubb """"7" w (8) 



f T _ TT . A cos Sab ^PaoPbo 1 . . 

Uab = Uab H „ ■ (9J 

^ Pao Pbo 

A second matrix, denoted X, contains differential operators Qxi — these are 
the second-order differential operators obtained from linearizing the quantum 
potential: 

2m x V vPx J 2m x \ ^Pxo + epx t J 



K 



2 



2m x 



(Qxo(Pxo)+sQxi(pxi)). (11) 



The quantity Qxo(pxo) corresponds to the background value of the quantum 
pressure, and contributes only to the background equations of motion — it 
does not affect the fluctuations. Now in a general background 



l{PXl) = -Z< [ 3 2 V + V ) PXl, (12) 

2 I Pxo Pxo Pxo J 
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and we define the matrix X to be 

X = - 



Qai 

o s. 







(13) 



Given the background homogeneity that will be appropriate for later parts 
of the current discussion, this will ultimately simplify to 



Qxi(pXi) 



1 



2p 



A'O 



in which case 



X 



m.A PAO 





"IB Pbo 



(14) 



(15) 



Without transitions between the two hyperfine states, when A = 0, the 
matrix S only contains the coupling constants 5ij — ■» Uij/h. While E is inde- 
pendent of the energy of the perturbations, the importance of X increases with 
the energy of the perturbation. In the so-called hydrodynamic approximation 
X can be neglected, effectively X — * and S — ► S. 

Besides the interaction matrix, we also introduce a transition matrix, 



^_ 2A cosJab y/PiQ Pjo 



-1 -1 
-1 +1 



and a mass-density matrix, 
D = h 

The final step is to define two column vectors, 



> = 


" 


= h 


d A 


" 





Pbo 
m B - 










and 



— [8ai, 8biY 



P = [PAl,PBlY 



(16) 

(17) 

(18) 
(19) 



We then obtain two compact equations for the perturbation in the phases and 
densities: 



e = -s p-v-we + ee, 

p= - v • (Dve + pv) - a e - e T p . 



(20) 
(21) 



Here the background velocity matrix simply contains the two background 
velocities of each condensate, 



V 



VAO 

vbo 



(22) 
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with two possibly distinct background velocities, 

n 



VBQ 



niA 
H 

tub 



-V0AO, 



(23) 



-V0BO- 



Additionally we also introduce the matrix 0, which depends on the difference 
of the initial phases and is defined as 



e = 



A sin Sab 



(24) 



Now combine these two equations into one: 

d t {5- 1 8) = -d t (s' 1 V-Ve)- V(V a 



V 



D- V V) V0 
1 



A 



(25) 



K9+-{T a d a 9 + d a (T a e)} 



where the index a runs from 0-3 (that is, over both time and space), and we 
now define 



r* = s- 1 e-e T s- 1 , 
r = VS^Q - e T s- 1 v, 



(26) 
(27) 



and 
K 



(28) 



Note that the r a matrices are antisymmetric in field-space (A <-> B), while 
the matrix K is symmetric. Also, both r a — > and K — > as <5ab — > 0. 

Our first goal is to show that equation (25), which fundamentally describes 
quasi-particle excitations interacting with a condensed matter system in the 
mean-field approximation, can be given a physical and mathematical interpre- 
tation in terms of a classical background geometry for massless and massive 
particles propagating through an analogue spacetime [2, 3, 4, 19]. This anal- 
ogy only holds (at least in its cleanest form) in the so-called hydrodynamic 
limit which limit is directly correlated with the healing length which 

we shall now introduce. 



2.3 Healing length 

The differential operator Qxi that underlies the origin of the X contribution 
above is obtained by linearizing the quantum potential 
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Vq( P x) = - 



2m x 



Ipx 



(29) 



which appears in the Hamilton-Jacobi equation of the BEC flow. This quan- 
tum potential term is suppressed by the smallness of h, the comparative large- 
ness of mxi and for sufficiently uniform density profiles. But of course in any 
real system the density of a BEC must go to zero at the boundaries of its 
electro-magnetic trap (given that px = \ipx (x, t)\ 2 ). In a 1-component BEC 
the healing length characterizes the minimal distance over which the order 
parameter goes from zero to its bulk value. If the condensate density grows 
from zero to po within a distance £ the quantum potential term (non local) 
and the interaction energy (local) are respectively ^kinetic ~ /i 2 /(2to£ 2 ) and 
^-interaction ~ Airh 2 apo / m. These two terms are comparable when 



£= (STrpoa)- 1 / 2 , 



(30) 



where a is the s-wave scattering length defined as 

m Un 



Airfr 2 



(31) 



Note that what we call Uq in the above expression is just the coefficient of the 
non- linear self-coupling term in the Gross-Pitaevskii equation, i.e., just L^aa 
or Ubb if we completely decouple the 2 BECs (Uab = A = 0). 

Only for excitations with wavelengths much larger than the healing length 
is the effect of the quantum potential negligible. This is called the hydrody- 
namic limit because the single-BEC dynamics is then described by the con- 
tinuity and Hamilton-Jacobi equations of a super-fluid, and its excitations 
behave like massless phonemic modes. In the case of excitations with wave- 
lengths comparable with the healing length this approximation is no longer 
appropriate and deviations from phononic behaviour will arise. 

Such a simple discrimination between different regimes is lost once one 
considers a system formed by two coupled Bosc-Einstcin condensates. One 
is forced to introduce a generalization of the healing £ length in the form 
of a "healing matrix" . If we apply the same reasoning used above for the 
definition of the "healing length" to the 2-componcnt BEC system we again 
find a functional form like that of equation (30) however we now have the 
crucial difference that both the density and the scattering length are replaced 
by matrices. In particular, we generalize the scattering length a to the matrix 
A: 



A 



1 



Uaa Uab 
Uab Ubb 



(32) 



Furthermore, from (30) a healing length matrix Y can be defined by 



Y~ 2 = — 

h 2 



y/PAorriA 
V^bo^b 



Uak Uab 
Uab Ubb 



V/°ao™a 

V*°B0 m B 



(33) 
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That is, in terms of the matrices we have so far defined: 

y-2 = 1 r i/2 SI -i/2. Y 2 = 2X 1 ' 2 ~- 1 X 1 ' 2 . (34) 

Define "effective" scattering lengths and healing lengths for the 2-BEC system 
as 

1 , r .i m A U AA + itlbUbb fo -s 

° cff = 2 ^ = 8^ ' (35) 



and 



f2 _ 1 . rv 2i _ , rv—n - h [^bb/(w a Pao) + ^aa/(^bPbo)] 
That is 

^2 _ fr 2 [mAPAoUAA + m B pBoUBB] 

4m A m-B/OA0/OB0 (^aa%b - U AB ) 

Note that if the two components are decoupled and tuned to be equivalent to 
each other, then these effective scattering and healing lengths reduce to the 
standard one-component results. 



3 Emergent spacetime at low energies 

The basic idea behind analogue models is to re-cast the equation for excita- 
tions in a fluid into the equation describing a masslcss or massive scalar field 
embedded in a pseudo— Riemannian geometry. Starting from a two component 
supcrfluid we are going to show that it is not only possible to obtain a massive 
scalar field from such an analogue model, in addition we are also able to model 
much more complex geometries. In Fig. 2 we illustrate how excitations in a 
2-component BEC are associated with various types of emergent geometry. 

Most generally, we show that excitations in a 2-component BEC (in the 
hydrodynamic limit) can be viewed as propagating through a specific class of 
pseudo-Finsler geometry. As additional constraints are placed on the BEC pa- 
rameter space, the geometry changes from pseudo-Finsler , first to bi-metric, 
and finally to mono-metric (pseudo Riemannian, Lorcntzian) geometry. This 
can be accomplished by tuning the various BEC parameters, such as the tran- 
sition rate A, the background velocities vao, i'boj the background densities 
Pao-, /°B0: and the coupling between the atoms Uaa, Ubb and Uab- 

At first, it might seem to be quite an artificial thing to impose such con- 
straints onto the system. But if one considers that the two macroscopic wave 
functions represent two interacting classical fields, it is more or less obvious 
that this is the only way in which to enforce physical constraints onto the 
fields themselves, and on the way they communicate with each other. 
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fab — rfTj ba 

pseudo-Finsler 




pad _ [j-T] ai > 


bi-metric ' 


Uab = 

_T 

(7ab = 
'11 




1 n 

QJ 

1 S "AO = «B0 

1 jjl d A = d B 
1 g Uab = 
d A C/ AA = d B U BB | g Uaa — Ubb 
1 

mono- metric | ' 







Fig. 2. The dependence of the emergent geometry on the 2-component BEC pa- 
rameters. 



3.1 Pseudo-Finsler geometry 

In the hydrodynamic limit (5 — > £), it is possible to simplify equation 
(25) — without enforcing any constraints on the BEC parameters — if we 
adopt a (3+l)-dimensional "spacetime" notation by writing x a = (t,x l ), with 
i € {1,2,3} and a £ {0,1,2,3}. Then equation (25) can be very compactly 
rewritten as [2, 3]: 

d a (f ab d b e) + (A + k) e + i {r a d a e + d a (v a e)} = o. (38) 

The object i ab is a 4 x 4 spacetime matrix (actually a tensor density), each of 
whose components is a 2 x 2 matrix in field-space — equivalently this can be 
viewed as a 2 x 2 matrix in field-space each of whose components is a 4 x 4 
spacetime tensor density. By inspection this is a self-adjoint second-order lin- 
ear system of PDEs. The spacetime geometry is encoded in the leading-symbol 
of the PDEs, namely the f ab , without considering the other subdominant 
terms. That this is a sensible point of view is most easily seem by considering 
the usual curved-spacetime d'Alembertian equation for a charged particle in- 
teracting with a scalar potential in a standard pseudo-Riemannian geometry 

-L[d a - iA a ] ( y/=g~g ab [d b - iA b ]e) +V9 = Q (39) 



from which it is clear that we want to make the analogy 
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f ab ~ ^—g g ab 
as the key quantity specifying the geometry. In addition 

r a ~ iA a and A + K ~ V - g ab A a A b 



(40) 



(41) 



so that T a is analogous to a vector potential and A (plus corrections) is related 
to the scalar potential V — in a translation invariant background this will 
ultimately provide a mass term. 

Specifically in the current 2-BEC system we have 





-(v~-T 


-vs- 1 


D-VS~ 1 V T 



where 



J A0 u 

vl 



(42) 



(43) 



is a 2 x 2 matrix in field space that is also a row vector in physical 3-spacc. 
Overall, this does look like a rather complicated object. However, it is possible 
to re-write the 4x4 geometry containing 2x2 matrices as its elements, in 
form of a single (2 • 4) x (2 • 4) matrix. 6 Explicitly we have 7 



pab 



: 11 



1 



-VA0 



—V 



AO 



"21 

which we can re- write as 
f 

where 





V A0 


VBO 





7 -l 

J 12 





«B0 


VAO 





"22 _ 



— V 



BO 



- «B0«B0 



(44) 





+ D xl h 


-5$ ViV 2 T - 


\-D 12 h~ 




+ D 21 h 


_ "22 i "^2V 2 i - 


\-D 22 h 



and 



■a b 



V? := (1,< ), 
V 2 a := (M B o), 

diag(0, 1,1,1) 



(45) 



(46) 
(47) 

(48) 



This result can be generalized for n-component systems. Any 4x4 geometry 
obtained from a n-component system can be re-written as a single {n ■ 4) X {n ■ 4) 
matrix. 

7 Note 5x2 = £21, so E~2 = . 
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Even simpler is the form 



f = 




"IS 








The key point is that this allows us 


to write 








fab 
ill 


-fab " 
/l2 








fab 
J21 


fab 
J 22 _ 


i 



where 



D®h. 



rab 

Jll — 


-7-1 
"11 




-D n h ab 


rab 

Jl2 — 


— -1 
"12 






rab 

J 21 — 


— -1 
"12 






rab 

J 22 — 


c-i 

"22 


V 2 a V 2 b 4 


- D 22 h ab 



(49) 



(50) 



(51) 



It is also possible to separate the representation of f ab into field space and 
position space as follows 



-n 1 





s 



-1 



Sf/ 




22 



Sal 1 



(52) 



V 2 Q V? 



a/) 



Why do we assert that the quantity f ab defines a pseudo-Finsler geometry? 
(Rather than, say, simply a 2 x 2 matrix of ordinary Lorentzian geometries?) 
To see the reason for this claim, recall the standard result [20] that the leading 
symbol of a system of PDEs determines the "signal speed" (equivalently, the 
characteristics, or the causal structure) [3]. Indeed if we consider the eikonal 
approximation (while still remaining in the realm of validity of the hydrody- 
namic approximation) then the causal structure is completely determined by 
the leading term in the Frcsncl equation 



dct[f ab fc a fc b ] = 0, 



(53) 



where the determinant is taken in field space. (The quantity f ab /c a fcb is exactly 
what is called the leading symbol of the system of PDEs, and the vanishing 
of this determinant is the statement that high-frequency modes can propa- 
gate with wave vector k a , thereby determining both characteristics and causal 
structure.) In the 2-BEC case we can explicitly expand the determinant con- 
dition as 

(fi b ik a kb){f^k c k d ) - (f?$k a k b )(fgk c k d ) = 0. (54) 
Define a completely symmetric rank four tensor 



Q 



abed 



r{ab red) 

111 J 22 



r(ab red) 
J 12 J 21 ' 



(55) 



Analogue spacetime based on 2-component Bose-Einstein condensates 



13 



then the determinant condition is equivalent to 

Q abcd k a k h k c k d = 0, (56) 

which now defines the characteristics in terms of the vanishing of the pscudo- 
co-Finsler structure 

Q(k) = Q abcd k a k b k c k c , (57) 

defined on the cotangent bundle. As explained in appendix A, this pscudo- 
co-Finsler structure can be Legendre transformed to provide a pseudo-Finslcr 
structure, a Finslerian notion of distance 

ds i = g a bcd dx a dx b dx c dx d . (58) 

Here the completely symmetric rank 4 tensor g a bcd determines the "sound 
cones" through the relation ds = 0. It is interesting to note that a distance 

function of the form 

ds = \J g a bcd dx a dx b dx c dx d (59) 

first made its appearance in Riemann's inaugural lecture of 1854 [21], though 
he did nothing further with it, leaving it to Finsler to develop the branch of 
geometry now bearing his name [22]. The present discussion is sufficient to 
justify the use of the term "pseudo-Finsler" in the generic 2-BEC situation, 
but we invite the more mathematically inclined reader to see appendix A for 
a sketch of how much further these ideas can be taken. 

The pseudo-Finsler geometry implicit in (50) is rather complicated com- 
pared with the pseudo-Riemannian geometry we actually appear to be living 
in, at least as long as one accepts standard general relativity as a good de- 
scription of reality. To mimic real gravity, we need to simplify our model. It 
is now time to use the major advantage of our analogue model, the ability to 
tune the BEC parameters, and with it the 2-ficld background configuration. 
The first order of business is to decouple f ab in field space. 

3.2 Bi-metric geometry 

The reduction of equation (52) to a diagonal representation in field space (via 
an orthogonal rotation on the fields), 

i ab - diag [/«* fg] = diag [V=gTi9il V=9^S&] , (60) 

enforces a bi-metric structure onto the condensate. There are two ways to 
proceed. 

Distinct background velocities 

For 



Vi + V 2 , 



(61) 
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we require all five 2x2 matrices appearing in (52) to commute with each 
other. This has the unique solution S^ 1 = 0, whence 

Uab = 0. (62) 

We then get 

"0 



~ n 1 o 





oh; 



V?V5 + Dh^. 



(63) 



Since D is a diagonal matrix this clearly represents a bi-mctric geometry. The 
relevant parameters arc summarized in Table 1. 



Equal background velocities 

For 

Vi = V 2 = V, (64) 

we are still dealing with a pseudo-Finsler geometry, one which is now indepen- 
dently symmetric in field space (f ab = [f T ] ab ), and position space i ab = f ba . 8 
In terms of the BEC parameters that means we must set equal the two back- 
ground velocities, v A q = v bo = v o, an d equation (52) is simplified to: 

f ab = _ s -l v a v b + m ab (65) 

From the above, diagonalizability in field space now additionally requires the 
commutator of the interaction and mass-density matrix to vanish: 

[3, D] = =>> U A B(d A -d B )=0. (66) 

Here, we have a choice between two tuning conditions that do the job: 

C/ AB = or d A =d B . (67) 

Under the first option, where Uab = 0, the two off-diagonal elements in 
equation (65) are simply zero, and we get the desired bi-metricity in the form 9 

V a V b + Dh ab . (68) 

Under the second option, for d A = d B = d, we have D = d I. The situation 
is now a bit trickier, in the sense that one has to diagonalize S : 



cab 



Sfi 1 o 

o e£ 



s- 1 = o T s~ 1 o 
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Bi-metric tuning scenarios 



VAO 7^ VBO 



vao = «bo 







Uab =0 



/it oc 



/22 OC 



"AO | AA A 



"0 | "AA a A ' 



„ „T 

"0"0 



u B0| (J BB a B ' 



"B0"B0 



"0 u BB d B ; 



--1 



gilab « 



S22 ai) OC 



-(c?i - vao) 


-vao 


-VAO 


ft" 


-(cf 2 - vlo) 




-«B0 


ft" 



-(cfi " 


vo 2 ) 




-V() 




h 13 


-(c 2 2 2 - 


vo 2 ) 




-v 


K 3 



C22 



Uaa<1a 



^AAPAO + ^ABPBO 
m A 



[>BBdB = ^BPBO + IWao 



Table 1. If the pseudo-Finsler geometry decouples into two independent Lorentzian 
geometries /ff — s/—gugu and /f* = y'— ffnffn, with two distinct speed of sounds 
en and C22, we are effectively dealing with a bi-metric Lorentzian metric. The table 
shows the results from three different tuning scenarios, that are sufficient to drive the 
2-component BEC from Finsler to bi-Lorentzian spacetime. The rightmost column 
dA = ds is addressed in [12] where the authors analyze cosmic inflation in such a 
bi-metric system. 



Once this is done, the way to proceed is to use the elements of S instead 
of EE -1 in equation (68). The relevant parameters are summarized in Table 1. 

There is a subtlety implicit in setting the background velocities equal that 
should be made explicit. If Vi = V-2 so that v A0 = ^bo, then since the masses 
appear in the relationship between phase and velocity we deduce 

m B A o(t, x) - m A 9 B o(t,x) = f(t). (70) 

If mA 7^ rn-B, and if the background velocity is nonzero, we must deduce 
that #ab(£, aO will be at the very least be position dependent, and we will 
be unable to set it to zero. Alternatively, if we demand <5ab = 0, and have 
V0ao(£, x ) = V0bo(£, x ) 7^ 0, then we cannot set uao = i>bo 7^ 0. Fortunately 
this will not seriously affect further developments. 



The most general pseudo-Finsler geometry is symmetric under simultaneous ex- 
change of field space and position space: f a = [f ] a . 

We would like to stress that this constraint can be easily fulfilled, at least in the 
special case Sab = 0, by tuning the transition rate A, see equation (9). 
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Last, but certainly not least, we present the conditions for a mono-metric 
geometry in a 2-component BEC. 



3.3 Mono-metric geometry 

Despite the fact that there are three different routes to bi-mctricity, once one 
demands mono-metricity, where 



f ab = diag [At, /#] = diag [V=Si7 Sir, Sit] 



(71) 



then one ends up with one set unique of constraints to reduce from pscudo- 
Finslcr to a single-metric Lorentzian geometry, namely: 



Uab = 0; 

Uaa = Ubb = U ; 

c?a = dji = d. 

This tuning completely specifies the spacctimc geometry, in that 

fab _ fab 
■I 11 — J22 cx 

and after a small calculation we get 



(72) 







-«o 


Udh l i - v v£, 



5ii 



ah 



922^ OC 





T 

-n 


-v 


h %3 



where we have defined 



Ud, 



(73) 



(74) 



(75) 



as the speed of sound. 10 

Throughout the preceding few pages we have analyzed in detail the first 
term in equation (38), and identified different condensate parameters with 
different emergent geometries. Since there is more then one term in the wave 
equation describing excitations in a two-component system, this is not the end 
of the story. The remaining terms in equation (38), which we might generically 
view as "mass" and "vector potential" terms, do not directly affect the space- 
time geometry as such. But when an excitation propagates through a specific 
analogue spacctime geometry, these terms will contribute to the kinematics. 
It then becomes useful to consider the "mass eigenmodes" in field-space. 



The speed of sound for quasi-particle excitations is of course our analogue for the 
speed of light in real gravity. 
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3.4 Merging spacetime geometry with mass eigenmodes 

The eigenmodes we are interested in are eigenmodes of the field-space matrices 
occurring in the sub-dominant terms of the wave equation. These eigenmodes 
(when they exist) do not notice the presence of multiple fields — in our specific 
case a 2-field system — and therefore propagate nicely through the effective 
curved spacetime. As promised in the abstract and the motivation, we are 
striving for an analogue model representing a massive scalar field in a mono- 
metric Lorentzian structure. By using the results from section 3.3 we are able 
to decouple the first term of equation (38). 

In the following we are focusing on two issues: First, we decouple the re- 
maining terms in equation (38), and subsequently we check that these eigen- 
modes do not recouple the geometric term. There is however one more (tech- 
nical) problem, and that is the fact that the terms we want to associate with 
the effective mass of the scalar field still contain partial derivatives in time 
and space, which ultimately implies a dependence on the energy of the prop- 
agating modes. 11 Luckily, this problem can be easily circumvented, for equal 
background phases, 12 

0ao = 0bo, (76) 

in which case 

K = T* = T 1 = 0. (77) 

This has the effect of retaining only the matrix A among the sub-dominant 
terms, so that the wave equation becomes 

d a ({ ab d b 6) + A8 = 0. (78) 

Due to the fact that the structure of the coupling matrix A cannot be changed, 
its eigenmodes determine the eigenmodes of the overall wave equation. The 
eigenvectors of A are given by 

EVl :=[+!, +1] 

EV2:=[-1,+1] v ; 

The final step is to make sure that our spacetime geometry commutes with 
the eigenvectors of A, that is 

[f afc ,vl]=0. (80) 

This constraint is only fulfilled in the mono-metric case, where we are dealing 
with two identical classical fields, that effectively do not communicate with 
each other. 13 That is, all field matrices are proportional to the identity matrix. 



11 This can be easily be seen by going to the eikonal approximation where d x — ► ik 
and dt — > iu. 

12 Note that <5ab = plus mono=metricity implies either mA = ibb with arbitrary 
vq 7^ 0, or tua 7^ rn-B with zero vq = 0. These are exactly the two situations we 
shall consider below. 

13 While Uab = 0, Uab + 0. 
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3.5 Special case: S — constant. 

There is one specific class of geometries we are particularly interested it, and 
that is when E is a position independent and time independent constant. In 
the next section we will focus exclusively on this case, and apply it to quantum 
gravity phenomenology. This case is however, also of interest as an example of 
an alternate interplay between fine tuning and emergent geometry. Under the 
assumption that E is position and time independent, we are able to directly 
manipulate the overall wave equation for the excitations and as a consequence 
obtain slightly milder tuning conditions for mono-metricity. 
Let us define 

6=~- 1/2 e, (8i) 

and multiply the whole wave equation (38) with E 1 ' 2 fr om the left. What we 
are doing is a transformation in field space onto a new basis 9, and in the new 
basis the wave equation is given by, 

d a (i ab d b 9) + (A + K) + i {f a dj + d a (f a 9)} = 0, (82) 

where the matrices in field space transform as: A = E x l 2 AS 1 / 2 , K = 
f a = E 1 / 2 T a E 1 / 2 , and the tensor-density as 



cab 



E l ' 2 i ab E x / 2 , (83) 



In general, the transformation matrix E 1 I 2 is a non-diagonal, though always 



symmetric: 14 



*V - (84) 
tr[H] + 2%/detS' 

A close look at equation (52), now using the tensor-density f ab , makes it 
obvious that for 

Uab = 0, (85) 

the geometry reduces from pseudo-Finsler to bi-metric. For the sake of keeping 
the discussion short and easy to follow, we set the background velocities equal, 
and now get 

fob = V ayb + p h ab (gg) 

In view of the tuning, Uab = 0, we see 

D = diag(J7 A A d A , Ubb d B ). (87) 

The new mass-density matrix, and therefore the overall geometry is diagonal 
in field space, hence we are now dealing with the required bi-metric structure. 



See appendix B. 
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So far we are in complete agreement with what we have obtained in our 
previous analysis, see Fig. 2. However, if we now ask for mono-metricity, we 
obtain a slightly milder constraint: 

Uaa d,A = Ubb cfe- 

Last but not least, we show in detail the results we obtain for this tuning 
scenario when including the A term (the mass term). To avoid confusion, we 
re-define a few matrices, 

C 2 = S 1 ' 2 D 5 1 / 2 ■ and Q 2 = S 1 ' 2 A S 1 ' 2 . (89) 

Both Cq and Q 2 are symmetric matrices. If [Cq, Q 2 ] — 0, which is equivalent 
to the matrix equation D E A — A E D, and is certainly satisfied in view 
of the above constraint, then they have common eigenvectors. Decomposition 
onto the eigenstates of the system results in a pair of independent Klein- 
Gordon equations 

* d a { v/=^ (gi /u ) ab d b e V n} + w?m 0i/ii = , (90) 

V -.9i/ii 1 > 

where the "acoustic metrics" arc given by 



(91/11)06 



v 



-v I I 



dxd 



(91) 



The metric components depend only on the background velocity Vq and the 
common speed of sound c. It is also possible to calculate the eigenfrequencies 
of the two phonon modes, 

cuj = 0; ujjj = tr[J7 2 ] . (92) 

A zero/ non-zero eigenfrequency corresponds to a zero/ non-zero mass for the 
phonon mode. 

In the eikonal limit we see that the in-phase perturbation will propagate 
with the speed of sound, 

v s = vq + k c , (93) 

while the anti-phase perturbations propagates with a lower group velocity 
given by: 

v g = ^r = v + k ==. 94 

ok yju jj + c 2 k 2 

Here k is the usual wave number. The dispersion relation we obtain for the 
mono- metric structure is Lorentz invariant. 

The fact that we have an analogue model representing both massive and 
massless particles is promising for quantum gravity phenomenology if we now 
extend the analysis to high-energy phonon modes where the quantum pres- 
sure term is significant, and where we consequently expect a breakdown of 
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Lorentz invariance. For the following, we concentrate on the generalization of 
flat Minkowski spacetimc, which implies a constant 5 and zero background 
velocities, v . In the language of condensed matter physics, we are thinking 
of a uniform condensate at rest. 



4 Application to quantum gravity phenomenology 

In using this 2-BEC model to probe issues of interest to the "quantum gravity 
phenomenology" community it behooves us to simplify as much as possible 
the parts of the model not of direct interest for current considerations. Specif- 
ically we wish to use the "quantum pressure" term as a model for the type 
of Lorentz violating physics that might occur in the physical universe at or 
near the Planck scale [23]. Since we are then interested in high energies, and 
consequently short distances, one might expect the average spacetime curva- 
ture to be negligible — that is, we will be interested in looking for "quantum 
pressure" induced deviations from special relativity, and can dispense with 
the notion of curved spacetimes for now. ("Flat" pseudo-Finsler spaces are 
already sufficiently complicated to lead to interesting physics.) In terms of the 
BEC condensates this means that in this section of the chapter we will con- 
centrate on a spatially-homogeneous time-independent background, so that 
in particular all the matrices f ab will be taken to be position- independent. 
(And similarly, S, A, D, etc. are taken to be position independent and we set 
vq = 0, so the background is at rest.) This greatly simplifies the calculations 
(though they are still relatively messy), but without sacrificing the essential 
pieces of the physics we are now interested in. 

Now the purpose of quantum gravity phenomenology is to analyze the 
physical consequences arising from various models of quantum gravity. One 
hope for obtaining an experimental grasp on quantum gravity is the generic 
prediction arising in many (but not all) quantum gravity models that ultra- 
violet physics at or near the Planck scale, M P i anck = 1.2 x 10 19 GeV/c 2 , 
(or in some models the string scale), typically induces violations of Lorentz 
invariance at lower scales [24, 25]. Interestingly most investigations, even if 
they arise from quite different fundamental physics, seem to converge on the 
prediction that the breakdown of Lorentz invariance can generically become 
manifest in the form of modified dispersion relations 



where the coefficients r\ n are dimensionless (and possibly dependent on the 
particle species considered), and we have restricted our expansion to CPT 
invariant terms (otherwise one would also get odd powers in k). The particular 
inertial frame for these dispersion relations is generally specified to be the 
frame set by cosmological microwave background, and ML or entz violation is the 




(95) 
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Fig. 3. How to tune the system to exhibit breakdown of Lorentz symmetry. 



scale of Lorentz symmetry breaking which furthermore is generally assumed 
to be of the order of Mpi anc k. 

Although several alternative scenarios have been considered in the liter- 
ature in order to justify the modified kinematics discussed above, to date 
the most commonly explored avenue is an effective field theory (EFT) ap- 
proach. In the present chapter we focus on the class of non-renormalizablc 
EFTs with Lorentz violations associated to dispersion relations like equa- 
tion (95). Relaxing our CPT invariance condition this class would include 
the model developed in [26], and subsequently studied by several authors, 
where an extension of quantum electrodynamics including only mass dimen- 
sion five Lorentz-violating operators was considered. (That ansatz leads to 
order k 3 Lorentz and CPT violating terms in the dispersion relation.) Very 
accurate constraints have been obtained for this model using a combination 
of experiments and observations (mainly in high energy astrophysics). See 
e.g. [25, 27, 28, 29]. In spite of the remarkable success of this framework as 
a "test theory" , it is interesting to note that there are still significant open 
issues concerning its theoretical foundations. Perhaps the most pressing one 
is the so called naturalness problem which can be expressed in the following 
way: Looking back at our ansatz (95) we can see that the lowest-order cor- 
rection, proportional to 772, is not explicitly Planck suppressed. This implies 
that such a term would always be dominant with respect to the higher-order 
ones and grossly incompatible with observations (given that we have very 
good constraints on the universality of the speed of light for different ele- 
mentary particles). Following the observational leads it has been therefore 
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often assumed either that some symmetry (other than Lorentz invariance) 
enforces the r) 2 coefficients to be exactly zero, or that the presence of some 
other characteristic EFT mass scale y, <C Mpi anc k (e.g., some particle physics 
mass scale) associated with the Lorentz symmetry breaking might enter in 
the lowest order dimcnsionlcss coefficient 772 — which will be then generically 
suppressed by appropriate ratios of this characteristic mass to the Planck 
mass: 772 oc (/Lt/Mpi anc k)' T where a > 1 is some positive power (often taken 
as one or two). If this is the case then one has two distinct regimes: For low 
momenta p/(Mpi anc kc) -C (/V-^Pianck)' 7 the lower-order (quadratic in the mo- 
mentum) deviations in (95) will dominate over the higher-order ones, while 
at high energies p/(Mpi anc kc) ^S> (/V^pianck)' 7 the higher order terms will be 
dominant. 

The naturalness problem arises because such a scenario is not well justified 
within an EFT framework; in other words there is no natural suppression of 
the low-order modifications in these models. In fact we implicitly assumed that 
there are no extra Planck suppressions hidden in the dimcnsionlcss coefficients 
rj n with n > 2. EFT cannot justify why only the dimcnsionlcss coefficients of 
the n < 2 terms should be suppressed by powers of the small ratio ///-Mpi anc k. 
Even worse, renormalization group arguments seem to imply that a similar 
mass ratio, /i/Aipi anc k would implicitly be present also in all the dimension- 
less n > 2 coefficients — hence suppressing them even further, to the point 
of complete undetectability. Furthermore it is easy to show [30] that, without 
some protecting symmetry, it is generic that radiative corrections due to par- 
ticle interactions in an EFT with only Lorentz violations of order n > 2 in 
(95) for the free particles, will generate n = 2 Lorentz violating terms in the 
dispersion relation, which will then be dominant. Observational evidence [24] 
suggests that for a variety of standard model particles |?/ 2 | ^ 10~ 21 . Natu- 
ralness in EFT would then imply that the higher order terms are at least as 
suppressed as this, and hence beyond observational reach. 

A second issue is that of "universality", which is not so much a "prob- 
lem", as an issue of debate as to the best strategy to adopt. In dealing with 
situations with multiple particles one has to choose between the case of univer- 
sal (particle-independent) Lorentz violating coefficients r/ n , or instead go for a 
more general ansatz and allow for particle-dependent coefficients; hence allow- 
ing different magnitudes of Lorentz symmetry violation for different particles 
even when considering the same order terms (same n) in the momentum ex- 
pansion. The two choices are equally represented in the extant literature (see 
e.g. [31] and [27] for the two alternative ansatze), but it would be interest- 
ing to understand how generic this universality might be, and what sort of 
processes might induce non-universal Lorentz violation for different particles. 

4.1 Specializing the wave equation 

For current purposes, where we wish to probe violations of Lorentz invariance 
in a flat analogue spacetime, we start with our basic wave equation (25) and 
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make the following specializations: i5ab — ¥ (so that r a — > and K — > 
0). We also set all background fields to be homogeneous (space and time 
independent), and use the formal operators 3 1 / 2 and 3~ x ^ 2 to define a new 
set of variables 

9 = S- 1 / 2 e, (96) 
in terms of which the wave equation becomes 

d 2 e = {3 1 / 2 [dv 2 -a\ s l ' 2 }e, (97) 

or more explicitly 

d 2 9 = { [3 - XV 2 ] 1 ' 2 [DV 2 -A] [3- XV 2 ] 1 / 2 } 9. (98) 

This is now a (relatively) simple PDE to analyze. The objects 3 1 ! 2 and 3~ x l 2 
are 2x2 matrices whose elements are pseudo-differential operators, but to 
simplify things it is computationally efficient to go directly to the eikonal 
limit where 15 

3^3 + Xk 2 . (99) 
This finally leads to a dispersion relation of the form 

det{w 2 1 - [3 + Xk 2 ] 1/2 [Dk 2 + A] [3 + Xk 2 ] 1 ' 2 } = , (100) 

and "all" we need to do for the purposes of this chapter, is to understand this 
quasiparticle excitation spectrum in detail. 



4.2 Hydro dynamic approximation 

The hydrodynamic limit consists of formally setting X — ► so that 3 — ► 3. 
(That is, one is formally setting the healing length matrix to zero: Y — > 0. 
More precisely, all components of the healing length matrix are assumed small 
compared to other length scales in the problem.) The wave equation (98) now 
takes the form: 

dfe = { - 1 / 2 [d\7 2 - A] 3 1 ' 2 } e. (ioi) 

Since this is second-order in both space and time derivatives, we now have at 
least the possibility of obtaining an exact "Lorentz invariance" . We can now 
define the matrices 

Q 2 = 3 1 ' 2 A 3 1 ' 2 - C 2 = 3 1 ' 2 D 3 1 ' 2 - (102) 



Once we are in the eikonal approximation the pseudo-differential operator 5 ' — > 
\/S + k 2 X can be given a simple and explicit meaning in terms of the Hamilton- 
Cay ley theorems of appendix B. 
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so that after Fourier transformation 

Lu 2 9 = {Cg k 2 + fl 2 } 6 = H{k 2 ) 6, (103) 
leading to the Fresnel equation 

det{w 2 I - H(k 2 )} = 0. (104) 

That is 

lu 4 - u? tr[H(k 2 )} + det[H(k 2 )} = 0, (105) 

whence 

j = tr[ff(fc 2 )] ± y/ggfgf - 4 dot[ff(F)T _ (iog) 

Note that the matrices f2 2 , C 2 , and H(k 2 ) have now carefully been arranged 
to be symmetric. This greatly simplifies the subsequent matrix algebra. Also 
note that the matrix H(k 2 ) is a function of k 2 ; this will forbid the appearance 
of odd powers of k in the dispersion relation — as should be expected due to 
the parity invariance of the system. 

Masses 

We read off the "masses" by looking at the special case of space-independent 
oscillations for which 

d 2 6 = -Q 2 0, (107) 

allowing us to identify the "mass" (more precisely, the natural oscillation 
frequency) as 

"masses" oc eigenvalues of (E 1 ! 2 A S 1 ^ 2 ) = eigenvalues of (5 A). (108) 

Since A is a singular 2x2 matrix this automatically implies 

w 2 = 0; J z II = tv{SA). (109) 

So we see that one mode will be a massless phonon while the other will have a 
non zero mass. Explicitly, in terms of the elements of the underlying matrices 

w 2 =0 ; u> 2 n = - VPA °/ B ° -{!7 AA + U BB - 2U AB ] (110) 



so that (before any fine-tuning or decoupling) 

AO / 
J? 



2 Pbo A 

= (HI) 



x J C/ AA + Ubb - 2U AB - 



y/PAO PBO 



Pao Pbo 

PBO V PAO 



It is easy to check that this quantity really does have the physical dimensions 
of a frequency 
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Mono-metricity conditions 

In order for our system to be a perfect analogue of special relativity: 

• we want each mode to have a quadratic dispersion relation; 

• we want each dispersion relation to have the same asymptotic slope. 

Let us start by noticing that the dispersion relation (106) is of the form 

lu 2 = [quadratic-,] ± \J [quartic] . (H2) 

The first condition implies that the quartic must be a perfect square 

[quartic] = [quadratic 2 ] 2 , (H3) 

but then the second condition implies that the slope of this quadratic must 
be zero. That is 

[quadratic 2 ](fc 2 ) = [quadratic 2 ](0), (H4) 

and so 

[quartic] (A: 2 ) = [quartic] (0) (115) 

must be constant independent of fc 2 , so that the dispersion relation is of the 
form 

uj 2 = [quadratic 1 ](fc 2 ) ± [quadratic 2 ](0). (H6) 

Note that this has the required form (two hyperbolae with the same asymp- 
totes, and possibly different intercepts). Now let us implement this directly in 
terms of the matrices Cq and M 2 . 

Step 1: Using the results of the appendix, specifically equation (255): 

dct[iJ 2 (fc)] = dct[^ 2 + C 2 k 2 } (117) 
= dct[n 2 } - tr {Q 2 C 2 } k 2 + dct[C 2 ] (fc 2 ) 2 . (118) 

(This holds for any linear combination of 2 x 2 matrices. Note that we apply 
trace reversal to the squared matrix Cg, we do not trace reverse and then 
square.) Since in particular det[,!? 2 ] = 0, we have: 

dct[ff 2 (fc)] = -tr {fl 2 Cl) k 2 + det[C 2 ] (fc 2 ) 2 . (119) 

Step 2: Now consider the discriminant (the quartic) 

quartic = tr[iJ(fc 2 )] 2 - 4 dct[iJ(fc 2 )] (120) 
= (tr[r? 2 ] + tr[C 2 ] fc 2 ) 2 - 4[ - tr {Q 2 C 2 } k 2 

+dct[C 2 ] (fc 2 ) 2 ] (121) 
= tr[i? 2 ] 2 + {2tr[X2 2 ]tr[C^] +4tr{/? 2 C$}}k 2 

+ |tr[C 2 ] 2 - 4dct[C 2 ]} (fc 2 ) 2 (122) 
= tr[J? 2 ] 2 + 2{2tr{Q 2 C$} - tr[/2 2 ]tr[C 2 ]}fc 2 

+ |tr[C 2 ] 2 - 4dct[C 2 ]} (fc 2 ) 2 . (123) 
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So in the end the two conditions above for mono-mctricity take the form 

( tr[C 2 ] 2 -4det[C 2 ] = 0: 
mono-metricity ^ j ^ fa ^ _ ^ ^ = Q (124) 

Once these two conditions are satisfied the dispersion relation is 
2 _ tr[H(k 2 )} ± tr[Q 2 ] _ tr[f2 2 ] ± tv[f2 2 ] + tr[C 2 ] k 2 



(125) 



whence 



1^ r^Oi ,9 9*9 9 r^9i 1 

2 



lu 2 = ^-tT[C£] k 2 = c 2 k 2 u'i = tr[f2 2 } + ^tr[C'£] k 2 = lu 2 h + c 2 k 2 , (126) 



as required. One mode is massless, one massive with exactly the "mass" pre- 
viously deduced. One can now define the quantity 

m n = Huju/c 2 ,, (127) 

which really does have the physical dimensions of a mass. 

Interpretation of the mono-metricity conditions 

But now we have to analyse the two simplification conditions 

CI : tr[C 2 ] 2 - 4 det[C* 2 ] = 0; (128) 
C2 : 2 tr {Q 2 C 2 } - tr[l2 2 ]tr[C 2 ] = 0; (129) 

to see what they tell us. The first of these conditions is equivalent to the 
statement that the 2x2 matrix Cq has two identical eigenvalues. But since Cq 
is symmetric this then implies Cq = c 2 I, in which case the second condition 
is automatically satisfied. (In contrast, condition C2 does not automatically 
imply condition CI.) Indeed if Cq = Cq I, then it is easy to see that (in order 
to make Cq diagonal) 

Uab = 0, (130) 



(which is sufficient, by itself, to imply bi-mctricity) and furthermore that 

Uaa Pao 2 U B b Pbo / 101 ^ 

= c o = ■ ( 131 ) 

m A m B 

Note that we can now solve for A to get 

A = -2^/pao Pbo Uab, (132) 

whence 

2 Uaa Pao + Uab Pbo U B b Pbo + Uab Pao on x 
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and 



J— + 


Ipso 


.V Pbo 


V PAO. 



*h = 4PA °; B 2 °^ AB |c/aa + U BB 2U AB + C/ AE 

(134) 

Note that (134) is equivalent to (112) with (132) enforced. But this then 
implies 



2 4 p aoPbqL/ab , 

J// = ^ 1 Uaa + ^ BB + ^ AB 



Paq Pbo 
Pbo Pao _ 



(135) 



Interpretation: Condition C2 forces the two low-momentum "propagation 
speeds" to be the same, that is, it forces the two 0(k 2 ) coefficients to be 
equal. Condition CI is the stronger statement that there is no 0(k 4 ) (or 
higher order) distortion to the relativistic dispersion relation. 

4.3 Beyond the hydro dynamical approximation 

At this point we want to consider the deviations from the previous analogue 
for special relativity. Our starting point is again equation (98), now retaining 
the quantum pressure term, which we Fourier transform to get: 

w 2 e = U/3 + Xk 2 [Dk 2 + A] y/S + Xk 2 \ 6 ' = H(k 2 ) 6. (136) 

This leads to the Fresnel equation 

det{w 2 I-H(k 2 )} = 0. (137) 



That is 
whence 



lu a - u? tr[H(k 2 )} + dct[H(k 2 )} = 0, (138) 



tr[H(k 2 )} ± y/ti[H{k 2 )] 2 ~4dct[H(k 2 )} 



u < = L v n v ^2 (139) 

which is now of the form 



u 2 = [quarticj ± y/ [octic] . (140) 



Masses 

The "masses" , defined as the zero momentum oscillation frequencies, are again 
easy to identify. Just note that the fc-independent term in the Fresnel equation 
is exactly the same mass matrix Q 2 = E 1 ! 2 A 5 1 / 2 that was present in 
the hydrodynamical limit. (That is, the quantum potential term X does not 
influence the masses.) 
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Dispersion relations 

Differently from the previous case, when the hydrodynamic approximation 
held, we now have that the discriminant of (139) generically can be an eighth- 
order polynomial in k. In this case we cannot hope to recover an exact analogue 
of special relativity, but instead can at best hope to obtain dispersion relations 
with vanishing or suppressed deviations from special relativity at low k; pos- 
sibly with large deviations from special relativity at high momenta. From the 
form of our equation it is clear that the Lorentz violation suppression should 
be somehow associated with the masses of the atoms m^m- Indeed we will use 
the underlying atomic masses to define our "Lorentz breaking scale" , which 
we shall then assume can be identified with the "quantum gravity scale" . The 
exact form and relative strengths of the higher-order terms will be controlled 
by tuning the 2-BEC system and will eventually decide the manifestation (or 
not) of the naturalness problem and of the universality issue. 

Our approach will again consist of considering derivatives of (139) in grow- 
ing even powers of k 2 (recall that odd powers of k are excluded by the parity 
invariance of the system) and then setting k — > 0. We shall compute only the 
coefficients up to order k as by simple dimensional arguments one can expect 
any higher order term will be further suppressed with respect to the A: 4 one. 

We can greatly simplify our calculations if before performing our analysis 
we rearrange our problem in the following way. First of all note that by the 
cyclic properties of trace 



tr[#(fc 2 )] = tr[(Dk 2 + A) (E + k 2 X)\ 

= tx[AS + k 2 {DE + AX) + [k 2 ) 2 DX] 
= tx[E x ' 2 AS 1 ' 2 + e^^DS 1 ' 2 + X 1 ' 2 AX 1 ' 2 ) 
+ {k 2 ) 2 X 1 l 2 DX 1 ' 2 ]. 



(141) 
(142) 



(143) 



Putting this all together, we can now define symmetric matrices 



Q 2 = E 1 / 2 AS 1 / 2 ] 

Cq = E 1 / 2 DE 1 ^ 2 ; AC 2 = X 1 ' 2 AX 1 ' 2 - 
C 2 = C 2 + AC 2 = E 1 ' 2 DE 1 ' 2 + X 1 ' 2 AX 1 ' 2 - 

h 2 

Z 2 = 2X 1/2 DX 1/2 = — M~ 2 . 

2 



(147) 



(144) 



(145) 
(146) 



With all these definitions we can then write 



tr[#(fc 2 )] = tr n 2 + k 2 {C 2 + AC 2 ) + -(k 2 ) 2 Z 2 



(148) 



where everything has been done inside the trace. If we now define 



H s (k 2 ) = Q 2 + k 2 (Cl + AC 2 ) + -(k 2 ) 2 Z 2 7 



(149) 
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then H s (k 2 ) is by definition both polynomial and symmetric and satisfies 

(150) 

while in contrast, 



tx[H(k 2 )]=tr[H s (k 2 )}, 
dct[H(k 2 )} ^ dct[H s (k 2 )}. 



to 



But then 



Whence 

dm 2 _ 1 
dfc 2 ~ 2 

and at k = 

dm 2 
dk 2 



tr[H s (k 2 )} ± yJtr[H s {k 2 )} 2 - 4det[H(k 2 )} 



tr[H'(k 2 )}± 



ti[H s (k 2 )]tr[H' s (k 2 )] - 2det'[ff(fc 2 )] 
y/tT[H s (k 2 )] 2 - 4dct[ff(F)J 



tr[C 2 ] ± 



tr[tt 2 ]tr[C 2 ] - 2dct'[H{k' 2 



Jfe^O 



But now let us consider 



tr[f2 2 



det[7J(fc 2 )] = det[(L>fc 2 + A) (S + k 2 X)} 
= dct[L>fc 2 + A] det[E + k 2 X] 



(151) 
(152) 

(153) 
(154) 



(155) 
(156) 



= det[sV 2 (Dk 2 + A)- 1 ' 2 ] det[J + k 2 S' 1 ' 2 XS' 1 ' 2 ] (157) 
where we have repeatedly used properties of the determinant. Furthermore 



dctLT + fc 2 ^ 1 / 2 ^- 1 / 2 



= det[I + k 2 S~ 1 X] 

= det[I + k 2 X 1/2 ~X 1/2 } 

= det[/ + fc 2 r 2 /2], 



so that we have 



det[H(k 2 )} = det[f2 2 + C$k 2 } det[7 + fc 2 F 2 /2] 



(158) 
(159) 
(160) 



(161) 



Note the the matrix Y 2 is the "healing length matrix" we had previously 
defined, and that the net result of this analysis is that the full determinant is 
the product of the determinant previously found in the hydrodynamic limit 
with a factor that depends on the product of wavenumber and healing length. 
But now, given our formula (255) for the determinant, we see 



det'[H(k 2 )} 



(-tr(^ 2 (7 2 ) + 2fc 2 det[C* 2 ]) det[J + fc 2 F 2 /2] 
+det[/? 2 + Clk 2 ] ("tr[? 2 ] + fc 2 dct[F 2 ])/2, (162) 



whence 
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dct'[£T(fc 2 )]^ = -tr(f2 2 C 2 ), 



and so 



That is: 



dm 2 
dk 2 

dm 2 
dk 2 



k^O 



k^O 



tr[C 2 } ± M^]tr[C 2 ] + 2tr(f2 2 C 2 ) 



tr[!2 2 ] 

t r[ c 2 ] ± {t r[ c 2 ] + 2 « 



(163) 
(164) 

(165) 



Note that all the relevant matrices have been carefully symmetrized. Also note 
the important distinction between C 2 and C 2 . Now define 



c 2 = ltr[C 2 ] 



then 



with 



dw 2 
dfc 2 



fc^O 

t2i i o±,.t cttni\ 



m 



f tr[C 2 ]tr[!2 2 ] +2tr(f2 2 C 2 ) 
{ tr[f2 2 ]tr[C 2 ] 
Similarly, consider the second derivative: 



= 1 + 



tr(l? 2 C 2 ) 
iv 2 n c 2 



(166) 
(167) 
(168) 



d(k 2 y 



t r [H'J(k 2 )} 

tr[H s (k 2 )}tr[H'J{k 2 )}+ti[H' s (k 2 )}tr[H' s {k 2 )} - 2det"[i?(fc 2 ) 



^/tr[H s (k 2 )} 2 -idct[H(k 2 )} 
(tr[H s {k 2 )}tr[H' s (k 2 )} - 2det' '[H {k 2 )}) 2 ' 



{tr[H s (k 2 )} 2 - 4det[i?(fc 2 )]) 3 / 2 



(169) 



whence 



d 2 ^ 2 




d(p) 2 


fe-»0 



2 tr[r? 2 ]tr[Z 2 ] + tr[C 2 ] 2 - 2det"[g(fc 2 )]^ 
U ' [Z 1 ± t^ 2 ] 

(tr[/2 2 ]tr[C 2 ]-2det'[i/(fc 2 )]^ ) 2 



tr[Sl 2 } 



(170) 



The last term above can be related to dw 2 /dfc 2 , while the determinant piece 
is evaluated using 

dct"[ff (fc 2 )] = (2det[C 2 ]) det[I + fc 2 F 2 /2] (171) 
+ (-tr(^ 2 C' 2 ) + 2fc 2 det[C 2 ]) (-tr[F 2 ] + fc 2 det[F 2 ])/2 
+dct[^ 2 + C 2 fc 2 ] (det[F 2 ]/2) 

+ (-tr(tt 2 C' 2 ) + 2fc 2 det[C 2 ]) (-tr[Y 2 ] + fc 2 det[F 2 ])/2. 
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Therefore 

det"[H{k 2 )) k ^o = (2dct[C 2 ]) 

+ (-tr(f2 2 C 2 )) (-tr[f 2 ])/2 + det[f2 2 ] (det[F 2 ])/2 

+ (-tr(^ 2 C 2 )) (-tr[F 2 ])/2. (172) 

That is, (recalling trL4] = -trL4]), 

det"[H(k 2 )] k ^o = (2dct[C 2 ]) - (tr(f2 2 C 2 )) (tr[F 2 ]), (173) 



or 



dct"[tf (fc 2 )] fc ^ = -tr[C 2 C- 2 ] - tr[/2 2 Cg] tr[y 2 
Now assembling all the pieces, a little algebra yields 

d 2 cc 2 



(174) 



d(fc 2 ) 2 



fe^O 



1 
1 



tr[^] ± tr[Z 2 ] ± 2^^gV 2 ] ± ^ " 



tr[!2 2 



tr[tt 2 



212 
21 



tr[/2 : 



(175) 



With the above formula we have completed our derivation of the lowest- 
order terms of the generic dispersion relation of a coupled 2-BEC system — 
including the terms introduced by the quantum potential at high wavenumber 
— up to terms of order k 4 . From the above formula it is clear that we do not 
generically have Lorentz invariance in this system: Lorentz violations arise 
both due to mode-mixing interactions (an effect which can persist in the 
hydrodynamic limit where Z — ► and Y — ► 0) and to the presence of the 
quantum potential (signaled by Z ^ and Y ^ 0). While the mode-mixing 
effects are relevant at all energies the latter effect characterizes the discrete 
structure of the effective spacetime at high energies. It is in this sense that 
the quantum potential determines the analogue of quantum gravity effects in 
our 2-BEC system. 



4.4 The relevance for quantum gravity phenomenology 

Following this physical insight we can now easily identify a regime that is 
potentially relevant for simulating the typical ansatze of quantum gravity 
phenomenology. We demand that any violation of Lorentz invariance present 
should be due to the microscopic structure of the effective spacetime. This 
implies that one has to tune the system in order to cancel exactly all those 
violations of Lorentz invariance which are solely due to mode-mixing interac- 
tions in the hydrodynamic limit. 

We basically follow the guiding idea that a good analogue of quantum- 
gravity-induced Lorentz violations should be characterized only by the ultra- 
violet physics of the effective spacetime. In the system at hand the ultraviolet 
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physics is indeed characterized by the quantum potential, whereas possible vi- 
olations of the Lorentz invariance in the hydrodynamical limit are low energy 
effects, even though they have their origin in the microscopic interactions. We 
therefore start by investigating the scenario in which the system is tuned in 
such a way that no violations of Lorentz invariance are present in the hydro- 
dynamic limit. This leads us to again enforce the conditions Cl and C2 which 
corresponded to "mono-metricity" in the hydrodynamic limit. 
In this case (165) and (175) take respectively the form 



dm 2 
dk 2 



= \ [tr[C 2 ] + (1 ± 1) tr[AC 2 }] = c 2 + I±l t r[zAC 2 ], (176) 



and 



d 2 L0 2 



t^l±t,^l Ttr[catt[r2 , 



1 tr[AC 2 } 2 + 2tr[C 2 ]tr[AC 2 } 1 tr[AC 2 ] 2 

2 trfft 5 ] T 2 tv[Q 2 ] 
tr[Z 2 ] ± tr[Z 2 ] ... / 2 tr[/AC 2 i 



d{k 2 



^^r^ + Wj' (17?) 

Recall (see section 4.2) that the first of the physical conditions Cl is equiv- 
alent to the statement that the 2x2 matrix C 2 has two identical eigenvalues. 
But since Cg is symmetric this then implies C 2 = Cg I, in which case the 
second condition is automatically satisfied. This also leads to the useful facts 

Uab = =► X = -2^ PA0 pbqUab; (178) 
2 £/aa Pao Ubb Pbo 



m A m B 



(179) 



Now that we have the fine tuning condition for the laser coupling we can 
compute the magnitude of the effective mass of the massive phonon and de- 
termine the values of the Lorentz violation coefficients. In particular we shall 
start checking that this regime allows for a real positive effective mass as 
needed for a suitable analogue model of quantum gravity phenomenology. 



Effective mass 

Remember that the definition of mjj reads 

mjj = h 2 L0 2 n /cQ. (180) 

Using equation (178) and equation (179) we can rewrite Cg in the following 
form 

eg = [mBPAo^AA + toaPbo^bb + £/ab(paowia + PboTOb)]/ (2to a tob)- (181) 
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Similarly equation (178) and equation (179) when inserted in equation (135) 
give 

2 _ 4[/ab(pao™b + 
Wj/ ^ . (182) 

We can now estimate m/j by simply inserting the above expressions in equa- 
tion (180) so that 

m 2 8L/ A b(pao"t-a + pbqWb )m A m B 

[m B pAo^AA + itiaPboUbb + C^ab(paowa + Pbo™b)] 

This formula is still a little clumsy but a great deal can be understood by 
doing the physically reasonable approximation m A ~ m-B = rn and p A ~ Pb- 
In fact in this case one obtains 



^ [U aa + TuZ + UbbY 



This formula now shows clearly that, as long as the mixing term [Tab is small 
compared to the "direct" scattering U aa + Ubb, the mass of the heavy phonon 
will be "small" compared to the mass of the atoms. Though experimental 
realizability of the system is not the primary focus of the current article, we 
point out that there is no obstruction in principle to tuning a 2-BEC system 
into a regime where \Uab\ <S! \U aa + Ubb\- For the purposes of this paper it 
is sufficient that a small effective phonon mass (small compared to the atomic 
masses which set the analogue quantum gravity scale) is obtainable for some 
arrangement of the microscopic parameters. We can now look separately at the 
coefficients of the quadratic and quartic Lorentz violations and then compare 
their relative strength in order to see if a situation like that envisaged by 
discussions of the naturalness problem is actually realized. 



Coefficient of the quadratic deviation 

One can easily see from (176) that the r/2 coefficients for this case take the 
form 

V2J = 0; (185) 
mji eg = tr[AC 2 ] = tT[X^ 2 AX^ 2 ] = tr[XA] 

1 A / m A p AQ + m B pBo 

2 m A mB \ VpaoPbo 

So if we insert the fine tuning condition for A, equation (178), we get 




C^ab (wapao + mB/OBo) 

mji = = 2 

niAniBCo 



(187) 



34 



Silke Weinfurtner, Stefano Liberati, and Matt Visser 



Remarkably we can now cast this coefficient in a much more suggestive form 
by expressing the coupling U AB m terms of the mass of the massive quasi- 
particle rrifj. In order to do this we start from equation (182) and note that 
it enables us to express U AB in (187) in terms ofwjj, thereby obtaining 



PAom-A + PBom B ujj 
772,7/ = 7-4 „ - : - ( 188 ) 



4c* p A0 m B + Pbo™a m A m B 

Now it is easy to see that 

Pao^a + PBom B 



Pao^b + Pbo™a 



0(1), (189) 



and that this factor is identically unity if cither m A — m B or pao — Pbo- All 
together wc arc left with 

m,u = v(^=) 2 , d90) 

V ^m A m B J 

where fj is a dimcnsionlcss coefficient of order unity. 

The product in the denominator of the above expression can be interpreted 
as the geometric mean of the fundamental bosons masses m A and m B . These 
are mass scales associated with the microphysics of the condensate — in anal- 
ogy with our experience with a 1-BEC system where the "quantum gravity 
scale" is set by the mass of the BEC atoms. It is then natural to define an ana- 
logue of the scale of the breakdown of Lorentz invariance as M e g = ^/m A m B . 
(Indeed this "analogue Lorentz breaking scale" will typically do double duty 
as an "analogue Planck mass".) 

Using this physical insight it should be clear that equation (190) effectively 
says 

which, given that mi — 0, we are naturally lead to generalize to 

/ mx \ 2 /mass scale of quasiparticle\ 2 _ ^ I II (192) 
\ M c ff J \ effective Planck scale / ' 

The above relation is exactly the sort of dimensionless ratio (p/M) CT that has 
been very often conjectured in the literature on quantum gravity phenomenol- 
ogy in order to explain the strong observational constraints on Lorentz viola- 
tions at the lowest orders. (See earlier discussion.) Does this now imply that 
this particular regime of our 2-BEC system will also show an analogue version 
of the naturalness problem? In order to answer this question we need to find 
the dimensionless coefficient for the quartic deviations, 774, and check if it will 
or won't itself be suppressed by some power of the small ratio mn/M e f[. 
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Let us now consider the coefficients of the quartic term presented in equa- 
tion (177). For the various terms appearing in (177) we get 



tr[Z 2 ] = 2tr[DX] = — 



m A + m\ 



2 9 
m A Tn B 



tr[AC 2 } = tr[XA] = -- 



A rriAPAo + niBPBo 



2 m A m B VPaoPbo 
tr[F 2 ] = 2tr[X~- 1 



= U A b 



niAP ao + m B p BQ 
m A m B 



h 2 PAom A UAA + PBam B U B B 



(193) 
(194) 
(195) 



2 PAomAPBomBUAAUBB 

where in the last expression we have used the fact that in the current scenario 
Uab = 0. Now by definition 



V4 = \(M 2 s /h 2 ) 



(dfc 2 



k=0 



is the dimcnsionless coefficient in front of the k . So 



Ml 



2h 2 



tr[Z 2 } ± tr[Z 2 ] 



±tr[C 2 ] 



M 2 r 2 
1V1 cS c 

h 2 



Whence 



tr[Z 2 } ± tr[Z 2 } 

MC 2 ] 

tr[Z 2 ] 



± 



M 2 r 2 
Ivl cS c 

H 2 
M 2 r 2 

h 2 



tr[C 2 ] 



tr[Y 2 } tr[AC 2 } 
2~ + tr[/2 2 ] 
tr[F 2 ] tr[Z\C 2 
2~ + tv[Q 2 ] 

tr[Y 2 ] tv[AC 2 
2 + tr[n 2 } 



tr[Y 2 } _ tr[AC 2 } 
2 tr[n 2 ] 



Let us compute the two relevant terms separately: 



tr[Z 2 ] 




f m\ + to 2 \ 




f m\ + m\ \ 


tr[C 2 ] 


4c 2 


V ™ A m B ) 


' 4c^M 2 ff 





(196) 

(197) 
(198) 

(199) 
(200) 

(201) 



-tr[F 2 ]/2 + 



tr [AC 2 
tr[f2 2 } 



AM 2 



PaoPbo^aa^bb [UaA + U BB 



AM 2 r 2 



m 2 A U aa + m B U BB 
mAmB (Uaa + Ubb 



(202) 
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where we have used pxoUxx = m xCo for X — A, B as in equation (179). 
Note that the quantity in square brackets in the last line is dimcnsionless. So 
in the end: 



mATOB 



mjjJAA + m B U BB 



m A m B 



(Uaa 



U 



BB 



m 2 A UBB + m^U AA 
m A m B [Uaa + U BB 



m 2 A U aa + "l|^BB 

rnAms (Uaa + U BB 



(203) 
(204) 
(205) 



Note: In the special case vtla = m B we recover identical quartic deviations 
7/4 j = 774.77 = 1/4, indicating in this special situation a "universal" deviation 
from Lorentz invariance. Indeed we also obtain 7/4 / = 7/4 /7 if we demand 
C^aa = U BB , even without fixing tua = m B . 

Thus in the analogue spacetime we have developed the issue of universal- 
ity is fundamentally related to the complexity of the underlying microscopic 
system. As long as we keep the two atomic masses ttia and m B distinct we 
generically have distinct 774 coefficients (and the 7/2 coefficients are unequal 
even in the case ttia = %)■ However we can easily recover identical 774 coef- 
ficients, for instance, as soon as we impose identical microphysics for the two 
BEC systems we couple. 



Avoidance of the naturalness problem 

We can now ask ourselves if there is, or is not, a naturalness problem present 
in our system. Are the dimensionless coefficients 7/4 7/77 suppressed below 
their naive values by some small ratio involving M e s = y/mAm B ? Or are 
these ratios unsuppressed? Indeed at first sight it might seem that further 
suppression is the case, since the square of the "effective Planck scale" seems to 
appear in the denominator of both the coefficients (204) and (205). However, 
the squares of the atomic masses also appear in the numerator, rendering both 
coefficients of order unity. 

It is perhaps easier to see this once the dependence of (204) and (205) on 
the effective coupling U is removed. We again use the substitution Uxx = 
it^xCq/ pxo for X — A, B, so obtaining: 



1 



mAPAo + raBPBO 
toaPbo + m B p ao 



m A p B0 



m B pAo 



m A m B (m A /JBo + "t-bPao) 



(206) 
(207) 
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From these expressions is clear that the 774,7/77 coefficients are actually of 
order unity. 

That is, if our system is set up so that 777,77 <C wa/b — which we have seen 
in this scenario is equivalent to requiring Uab *C £/aa/bb — no naturalness 
problem arises as for p > 777,77 cq the higher-order, energy-dependent Lorentz- 
violating terms (n > 4) will indeed dominate over the quadratic Lorentz- 
violating term. 

It is quite remarkable that the quadratic coefficients (192) are exactly of 
the form postulated in several works on non-renormalizable EFTwith Lorentz 
invariance violations (see e.g. [25]). They are indeed the squared ratio of the 
particle mass to the scale of Lorentz violation. Moreover we can see from 
(204) and (205) that there is no further suppression — after having pulled 
out a factor (fi/MLorentz violation) 2 — for the quartic coefficients 774,1/11- These 
coefficients arc of order one and generically non-universal, (though if desired 
they can be forced to be universal by additional and specific fine tuning). 

The suppression of 772, combined with the non-suppression of 774, is pre- 
cisely the statement that the "naturalness problem" does not arise in the 
current model. We stress this is not a "tree level" result as the dispersion re- 
lation was computed directly from the fundamental Hamiltonian and was not 
derived via any EFTreasoning. Moreover avoidance of the naturalness prob- 
lem is not directly related to the tuning of our system to reproduce special 
relativity in the hydrodynamic limit. In fact our conditions for recovering spe- 
cial relativity at low energies do not a priori fix the the 772 coefficient, as its 
strength after the "fine tuning" could still be large (even of order one) if the 
typical mass scale of the massive phonon is not well below the atomic mass 
scale. Instead the smallncss of 772 is directly related to the mass-generating 
mechanism. 

The key question is now: Why docs our model escape the naive predic- 
tions of dominant lowest-dimension Lorentz violations? (In fact in our model 
for any p 3> 777,77 the k A Lorentz violating term dominates over the order 
k 2 one.) We here propose a nice interpretation in terms of "emergent sym- 
metry" : Non-zero A simultaneously produces a non-zero mass for one of the 
phonons, and a corresponding non-zero Lorentz violation at order k 2 . (Single 
BEC systems have only fc 4 Lorentz violations as described by the Bogoliubov 
dispersion relation.) Let us now drive A — > 0, but keep the conditions CI and 
C2 valid at each stage. (This also requires C/ab — * 0.) One gets an EFT which 
at low energies describes two non-interacting phonons propagating on a com- 
mon background. (In fact 772—5-0 and C7 = C77 = cq.) This system possesses 
a SO(2) symmetry. Non-zero laser coupling A softly breaks this 50(2), the 
mass degeneracy, and low-energy Lorentz invariance. Such soft Lorentz viola- 
tion is then characterized (as usual in EFT) by the ratio of the scale of the 
symmetry breaking 777,77, and that of the scale originating the Lorentz viola- 
tion in first place Mi, OTentz violation- We stress that the 50(2) symmetry is an 
"emergent symmetry" as it is not preserved beyond the hydrodynamic limit: 
the 774 coefficients are in general different if ttia 7^ 7773, so 50(2) is generi- 
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cally broken at high energies. Nevertheless this is enough for the protection 
of the lowest-order Lorentz violating operators. The lesson to be drawn is 
that emergent symmetries are sufficient to minimize the amount of Lorentz 
violation in the lowest-dimension operators of the EFT. In this regard, it is 
intriguing to realise that an interpretation of SUSY as an accidental symmetry 
has indeed been considered in recent times [32], and that this is done at the 
cost of renouncing attempts to solve the hierarchy problem in the standard 
way. It might be that in this sense the smallness of the particle physics mass 
scales with respect to the Planck scale could be directly related to smallness 
of Lorentz violations in renormalizable operators of the low-energy effective 
field theory we live in. We hope to further investigate these issues in future 
work. 



5 Outlook, summary and discussion 

So where can (and should) we go from here? If 2-component BECs provide 
such a rich mathematical and physical structure, are 3-componcnt BECs, or 
general multi-component BECs even better? That depends on what you are 
trying to do: 

• If one wishes to actually build such an analogue spacetime in the labora- 
tory, and perform actual experiments, then iteration through 1-BEC and 
2-BEC systems seems the most promising route in terms of our technolog- 
ical capabilities. 

• For n-componcnt BECs we sketch the situation in figure 4. The key point 
is that due to overall translation invariance one again expects to find one 
massless quasi-particle, with now n — 1 distinct massive modes. Unfortu- 
nately the matrix algebra is now considerably messier — not intrinsically 
difficult (after all we are only dealing with n x n matrices in field space) — 
but extremely tedious. Physical insight remains largely intact, but (except 
in some specific particularly simple cases), computations rapidly become 
lost in a morass of technical detail. 

• However, if one wishes to draw general theoretical lessons from the ana- 
logue spacetime programme, then multi-component systems arc definitely 
the preferred route — though in this case it is probably better to be even 
more abstract, and to go beyond the specific details of BEC-based systems 
to deal with general hyperbolic systems of PDEs. 

• In appendix A we have sketched some of the key features of the pseudo- 
Finsler spacetimes that naturally emerge from considering the leading sym- 
bol of a hyperbolic system of PDEs. While it is clear that much more could 
be done based on this, and on extending the field theory "normal modes" 
of [2, 3], such an analysis would very much move outside the scope of the 
COSLAB programme. 
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In short the 2-BEC system is a good compromise between a system complex 
enough to exhibit a mass-generating mechanism, and still simple enough to 
be technologically tractable, with good prospects for laboratory realization of 
this system in the not too distant future. 
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Fig. 4. The figure captures the key features of possible eigenmodes for a small 
perturbation (circles) in a 1 (left side), 2, 3, 4, and 5-component (right side) BEC. 
In a 1-component system only one kind of perturbation is allowed, which corresponds 
to a massless particle propagating through an effective curved spacetime, while in 
the 2-component case two different kinds of mode appear, the one in-phase (massless 
particle) and one in anti-phase (massive particle). For a three-component system we 
again expect to find one massless particle, when all perturbations are in phase, and 
now in addition to that two massive particles. 



The key features we have emphasised in this chapter have been: 

• A general analysis of the 2-BEC system to see how perturbations on a 
2-BEC background lead to a system of coupled wave equations. 

• Extraction of the geometric notion of pseudo Finsler spacetime from this 
wave equation, coupled with an analysis of how to specialize pscudo- 
Finsler geometry first to a bi-metric Lorentzian geometry and finally to the 
usual mono-metric Lorentzian geometry of most direct interest in general 
relativity and cosmology. 

• The mass-generating mechanism we have identified in suitably coupled 2- 
componcnt BECs is an essential step in making this analogue spacetime 
more realistic; whatever one's views on the ultimate theory of "quantum 
gravity" , any realistic low-energy phenomenology must contain some mass- 
generating mechanism. 

• Use of the "quantum pressure" term in the 2-BEC system to mimic the 
sort of Lorcntz violating physics that (based on the relatively young field 
of "quantum gravity phenomenology") is widely expected to occur at or 
near the Planck scale. 

• Intriguingly, we have seen that in our specific model the mass-generating 
mechanism interacts with the Lorentz violating mechanism, naturally lead- 
ing to a situation where the Lorentz violations are suppressed by powers of 
the quasi-particle mass scale divided by the analogue of the Planck scale. 

In summary, while we do not personally believe that the real universe is 
an analogue spacetime, we arc certainly intrigued by the fact that so much of 
what is normally viewed as being specific to general relativity and/or particle 
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physics can be placed in this much wider context. We should also be forthright 
about the key weakness of analogue models as they currently stand: As we 
have seen, obtaining an analogue spacetimc geometry (including spacetimc 
curvature) is straightforward — but what is not straightforward is obtaining 
the Einstein equations. The analogue models are currently analogue models of 
quantum field theory on curved spacetime, but not (yet?) true analogue mod- 
els of Einstein gravity. Despite this limitation, what can be achieved through 
the analogue spacetimc programme is quite impressive, and we expect in- 
terest in this field, both theoretical and hopefully experimental, to continue 
unabated. 



A Finsler and co— Finsler geometries 

Finsler geometries are sufficiently unusual that a brief discussion is in order 
— especially in view of the fact that the needs of the physics community are 
often somewhat at odds with what the mathematical community might view 
as the most important issues. Below are some elementary results, where we 
emphasise that for the time being we are working with ordinary "Euclidean 
signature" Finsler geometry. For general references, see [33]. 

A.l Basics 

Euler theorem: If H(z) is homogeneous of degree n then 



dz 



nH{z). (208) 



Finsler function: Defined on the "slit tangent bundle" T^o(M) such that 
F : T^ (M) -► [0, +oo) where 

F(x,t): F(x,\t) = XF(x,t), (209) 

and 

TMM) = |J [T x - {0}] . (210) 

xeM 

That is, the Finsler function is a defined only for nonzero tangent vectors 
t G [T x — {0}], and for any fixed direction is linear in the size of the vector. 
Finsler distance: 



<Ly(x,y) 
Finsler metric: 



rv 

I F (x{t) , dx / dr) dr; r= arbitrary parameter. (211) 

J X 



^ (M) = 2 WW ■ (212) 
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The first slightly unusual item is the introduction of co-Finsler structure: 
co-Finsler function: Define a co-Finsler structure on the cotangent bundle by 
Legendre transformation of F 2 (x,t). That is: 



G 2 (x,p)=t j (p) Pj -F 2 (x,t(p)) 
where t(p) is defined by the Legendre transformation condition 



d[F 2 
dP 



(x,t) =pj. 



Note 



dPj 
dt k 



d[F 2 



dP dt k 

which is why we demand the Finsler metric be nonsingular. 
Lemma: G(x,p) defined in this way is homogeneous of degree 1. 
Proof: Note 

, dH(z) 



(213) 



(214) 



(215) 



dz 



implies 



d_ 



(dz)> 



;H(Z) 



nH{z) 



= (n — in) 



(d z y 



Mz) 



(216) 



(217) 



In particular: 

F 2 is homogeneous of degree 2. 
gij is homogeneous of degree 0. 
d[F 2 ]/dt is homogeneous of degree 1. 
Therefore p(t) is homogeneous of degree 1 
and t{p) is homogeneous of degree 1. 
Therefore t{p)p — F 2 (t{p)) is homogeneous of degree 2. 
Therefore G(p) is homogeneous of degree 1. 

Thus from a Finsler function F(x, t) we can always construct a co-Finsler 
function G(x,p) which is homogeneous of degree 1 on the cotangent bundle. 

From the way the proof is set up it is clearly reversible — if you are given 
a co-Finsler function G{x,p) on the cotangent bundle this provides a natural 
way of extracting the corresponding Finsler function: 



F 2 (x,t) = tp(t)-G 2 (x,p(t)). 



(218) 



A. 2 Connection with the quasi-particle PDE analysis 

From the PDE-bascd analysis we obtain the second-order system of PDEs 

d a (f ab AB d b 9 B ) + lower order terms = 0. (219) 

We are now generalizing in the obvious manner to any arbitrary number n 
of interacting BECs. but the analysis is even more general than that — it 
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applies to any field-theory normal-mode analysis that arises from a wide class 
of Lagrangian based systems [2, 3]. 

Going to the eikonal approximation this becomes 

/ab PaPb e B + lower-order terms = 0, (220) 

which leads (neglecting lower order terms for now) to the Fresnel-like equation 

det[/f B p aPb ] = 0. (221) 

But by expanding the n x n determinant (n is the number of fields, not the 
dimension of spacctime) we have 

det[/X h B p aPb ] = Q abcd - PaPbPcPd ■ ■ ■ (222) 

where if there are n fields there will be 2n factors of p. 
Now define 

Q(x, P ) =Q abcd - PaPbPcPd ..., (223) 

and 

G(x,p) = 2 ^QM) = [Q{z,p)] 1/(2n \ (224) 

then 

• Q(x,p) is homogeneous of degree 2n. 

• G(x,p) is homogeneous of degree 1, and hence is a co-Finsler function. 

• We can now Legendre transform G — > F, providing a chain 

Q(x,p)^G(x,p)^F(x,t). (225) 

Can this route be reversed? 

Step 1: We can always reverse F(x,t) — > G{x,p) by Legendre transformation. 
Step 2: We can always define 

^-\wA [G{x ^ pn (226) 

this is homogeneous of degree 0, but is gencrically not smooth at p = 0. 
In fact, if g ab (x,p) is smooth at p = then there exits a limit 

g ab (x lP ^0)=g a \x), (227) 

but since g ab (x,p) is homogeneous of degree this implies 

g ab (x,p)=g ab (x) [Vp], (228) 

and so the geometry simplifies Finsler — » Riemann. 

This observation suggests the following definition. 
Definition: A co-Finsler function G(x, P ) is 2n-smooth iff the limit 
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(2n)iaUJ <%*>*r = Cr~~ (229) 

exists independent of the direction p in which you approach zero. 
Lemma: If G(x,p) is 2n-smooth then 

G(x, P ) 2n = Q abcd - PaPbPcPd . . . , (230) 

and indeed 

G(X, P ) = 2 \lQ abcd - PaPbPcPd--- (231) 

Proof: G 2n is homogeneous of degree 2n, so (d/dp) 2n G 2n is homogeneous of 
degree 0. Therefore if the limit 

/ r\ \ In 
1 / O 



(2n)! p^o 
exists, it follows that 



lim — G{x,pY n = Q abcd - (232) 



1 / O 



(2n)! 



G(x,p) 211 = Q o6cd - [Vp], (233) 



and so the result follows. 

Special case n = 1; If G(x,p) is 2-smooth then 

^^k G( ^ rf2 -^-^ p) ' (234) 

and co-Finslcr — > Riemann. 

These observations have a number of implications: 

• For all those co-Finsler functions that are 2n smooth we can recover the 
tensor Q abcd -. 

• Not all co-Finsler functions are 2n smooth, and for those functions we 
cannot extract Q abcd - - hi any meaningful way. 

• But those specific co-Finsler functions that arise from the leading symbol 
of a 2nd-order system of PDEs are naturally 2n-smooth, and so for the 
specific co-Finsler structures we are physically interested in 

Q{x,p) <-> G(x,p) <-> F(x, t). (235) 

• Therefore, in the physically interesting case the Finsler function F(x, t) 
encodes all the information present in Q abcd - - . 

Special case n = 2: For two fields (appropriate for our 2-BEC system), we can 
follow the chain 

f ab -> Q(x,p) <-> G(x,p) <-> F(x, t) (236) 

to formally write 
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ds i 



9abcd 



dx q dx b dx c dx' 



(237) 



ol- 



ds = \/ 1 g a bcd dx q dx b dx c dx d . 



(238) 



This is one of the "more general" cases Riemann alludes to in his inaugural 
lecture of 1854 [21]. 

This discussion makes it clear that the general geometry in our 2-BEC 
system is a 4-smooth Finsler geometry. It is only for certain special cases that 
the Finsler geometry specializes first to "multi-metric" and then to "mono- 
metric" Riemannian geometries. 

A. 3 Lorentzian signature Finsler geometries 

The distinction between Finsler and pseudo-Finsler geometries has to do with 
the distinction between elliptic and hyperbolic PDEs. Elliptic PDEs lead to 
ordinary Finsler geometries, hyperbolic PDEs lead to pseudo-Finsler geome- 
tries. 

Remember that in special relativity we typically define 



• d-y(x,y) G M + for spacelike paths; 

• d-y(x, y) — for null paths; 

• d~ l {x,y) G S + for timclike paths; 

The point is that even in special relativity (and by implication in general 
relativity) "distances" do not have to be real numbers. This is why physicists 
deal with pseMdo-Ricmannian [Lorentzian] geometries, not (strictly speaking) 
Riemannian geometries. 

To see how this generalizes in a Finsler situation let us first consider a 
co-Finsler structure that is multi-metric, that is: 



where each one of these n factors contains a Lorentzian signature matrix and 
so can pass through zero. Then 




(239) 



then 



Q(x,p) = n? =1 (g?>p a p b ), 



(240) 




(241) 



and 




(242) 



where 
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• I = — > G(x,p) £ M + — > outside all n signal cones; 

• ^ = n^G(a;,p)e J + — > inside all n signal cones. 

So we can now define 

• Spacelike «-» outside all n signal cones <-> G real; 

• Null <-> on any one of the n signal cones <-> G zero; 

• Timelike <-> inside all n signal cones <-> G imaginary; 

• plus the various "intermediate" cases: 

"intermediate" «-> inside £ oi n signal cones <-> G G / /n x 1? 



(243) 



Now this basic idea survives even if we do not have a multi-metric theory. 
The condition Q(x,p) = defines a polynomial of degree 2n, and so de- 
fines n nested sheets (possibly crossing in places). Compare with Courant and 
Hilbert's discussion of the Monge cone [20]. 
That is: 

Q(x,p)=0^Q(x,(E,p)) = 0; 

<^> polynomial of degree 2n in E for any fixed p; 

in each direction 3 2n roots in E; 
O corresponds to n [topological] cones. 

(These are topological cones, not geometrical cones, and the roots might hap- 
pen to be degenerate.) 

Question: Should we be worried by the fact that the co-metric g ab is singular 
on the signal cone? (In fact on all n of the signal cones.) Not really. We have 



G(x,p) = VQ 



~)abcd. . 



PaPbPcPd 



(244) 



so 



9 ab (x,p) 



whence 



1 d 2 



2 dp a dp b 



( \/Q{x,p) 



1 d 

2n dpb 



{Q^- 1 Q abcd -p bPcPd ...}, 

(245) 



1 



g a \x,p) = — Q™ Q a 



bed. . 



i /i 



2n 

1 

2n \n 

which we can write as 



PcPd 



i)Q^- 2 [Q 



^acde... 



PcPdPe 



9 



ab (x, P ) = — g^^"- 1 )/" Q abcd - 



2n 



PcPd 



1 n- 1 
2n n 



Q 



-(2n-l)/r. 



[Q 



a cde 



PcPdPe 



(246) 



[Q 



bfgh... 



PfPgPh ■ 



(247) 



[Q 



bfgh. 



PfPgPh ■ ■ 
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Yes, this naively looks like it's singular on the signal cone where Q(x,p) = 0. 
But no, this is not a problem: Consider 

9 ab Pa Pb = ^Q- (2n - 1)/n Q ^—Q- {2n - 1)/n Q 2 , (248) 

then 

and this quantity is definitely non-singular. 

A. 4 Summary 

In short: 

• pseudo-Finsler functions arise naturally from the leading symbol of hy- 
perbolic systems of PDEs; 

• pseudo-Finsler geometries provide the natural "geometric" interpretation 
of a multi-component PDE before fine tuning; 

• In particular the natural geometric interpretation of our 2-BEC model 
(before fine tuning) is as a 4-smooth pseudo-Finsler geometry. 

B Some matrix identities 

To simplify the flow of argument in the body of the paper, here we collect a 
few basic results on 2 x 2 matrices that are used in our analysis. 

B. l Determinants 

Theorem: For any two 2x2 matrix A: 

det(A) = i {tr[A} 2 - tr[A 2 } } . (250) 

This is best proved by simply noting 

dct(A) = AxA 2 = l - [(Ax + A 2 ) 2 - (A? + A 2 )] = \ {trL4] 2 - trL4 2 ]} . (251) 

If we now define 2x2 "trace reversal" (in a manner reminiscent of standard 
GR) by 

A = A - trL4] I; A = A; (252) 
then this looks even simpler 

det(A) = -^tr[A A] = det(A). (253) 
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A simple implication is now: 

Theorem: For any two 2x2 matrices A and B: 

det(A + A B) = det(A) + A {tr[A]tr[B] - tr[A B}} + A 2 dct(B). (254) 

which we can also write as 

det(A + A B) = det(A) - A tx[A B] + A 2 dct(B). (255) 

Note that tr[A B] = tx[A B\. 

B.2 Hamilton Cayley theorems 

Theorem: For any two 2x2 matrix A: 

tv[A] I A A_ 

A ~ det[A] " det[A]' (255j 

Theorem: For any two 2x2 matrix A: 



^- I \ A± ^ 1 1 (257 ) 
'tr[A]±2\/detI J 
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